Relativistic Kinetic Equations for Finite Domains and Freeze-out Problem 
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The relativistic kinetic equations for the two domains separated by the hypersurface with both 
space- and time-like parts are derived. The particle exchange between the domains separated by 
the time-like boundaries generates source terms and modifies the collision term of the kinetic equa- 
tion. The correct hydrodynamic equations for the "hydro+cascade" models are obtained and their 
differences from existing freeze-out models of the hadronic matter are discussed. 
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I. INTRODUCTION 

In recent years an essential progress has been achieved 
in our understanding of the freeze-out problem in rel- 
ativistic hydrodynamics, i.e. how to convert hydrody- 
namic solution into free streaming particles. Thus, in 
works Q the correct generalization of the famous Cooper- 
Frey formula S for the time-like hypersurfaces, the cut- 
off formula, was derived for the first time. This was a 
necessary, but yet not sufficient step to formulate the rel- 
ativistic hydrodynamics without causal paradoxes. The 
main problem was to formulate the energy-momentum 
and particle conservation not only for the expanding fluid 
alone, but to extend it onto the system consisting of the 
fluid and the gas of free streaming particles which are 
emitted from the freeze-out hypersurface. The princi- 
pal solution of this task was given in 0] and its further 
analysis is presented in ||. The advantage of this ap- 
proach is the absence of the logical and causal paradoxes 
usually arising, if emission of particles happens from the 
time-like hypersurfaces The disadvantage of this 

approach lies in its tremendous numerical complexity. 

Numerous attempts to improve and to develop this ap- 
proach further using primitive kinetic models were 
not very successful so far. Very recently a more funda- 
mental treatment || based on the analysis of the Boltz- 
mann equation was suggested. This approach, however, 
did not overcome the usual difficulty of transport equa- 
tions in describing the phase transition phenomenon. 
This difficulty has been overcome naturally within the 
"hydro + cascade" models suggested in Ref. JtJ (BD) 
and further developed in || (TLS). The latter mod- 
els assume that the nucleus-nucleus collisions proceed 
in three stages: hydrodynamic expansion ("hydro") of 
quark gluon plasma (QGP), phase transition from QGP 
to hadron gas (HG) and the stage of hadronic rescatter- 
ing and resonance decays ("cascade"). The switch from 
hydro to cascade modeling takes place at the boundary 
between the mixed and hadronic phases. The spectrum of 
hadrons leaving this hypersurface of the QGP-HG tran- 
sition is taken as input for the cascade. 

Evidently, such an approach incorporates the most at- 
tractive features of both hydrodynamics, which describes 
the QGP-HG phase transition very well, and cascade, 



which works better during hadronic rescattering. How- 
ever, both the BD and TLS models face some principal 
difficulties which cannot be ignored. Thus, within the 
BD approach the initial distribution for cascade is found 
by the Cooper- Frye formula , which takes into account 
particles with all possible velocities, whereas in the TLS 
model the initial cascade distribution is given by the cut- 
off formula jlB , which accounts for only those particles 
that can leave the phase boundary. As shown below the 
Cooper-Frye formula will lead to causal and mathemati- 
cal problems in the present version of BD model because 
the QGP-HG phase boundary inevitably has time-like 
parts. On the other hand the TLS model from the be- 
ginning does not conserve energy, momentum and num- 
ber of charges and this, as demonstrated later ,is due to 
the fact that the equations of motion used in [|8| are not 
complete and, hence, should be modified. 

The main difficulty of the "hydro + cascade" approach 
looks very similar to the freeze-out problem in relativistic 
hydrodynamics. In both cases the finite domains (sub- 
systems) have time- like boundaries through which the ex- 
change of particles is occurring and this should be taken 
into account. In relativistic hydrodynamics the prob- 
lem was solved by the constraints which appeared on the 
freeze-out hypersurface and provided the global energy- 
momentum and charge conservation Similarly, in 
kinetic theory one has to modify the transport equa- 
tions by the source terms which describe the exchange 
of the particles on the time-like parts of the boundary 
between domains. Therefore, we shall consider the two 
semi-infinite domains which are separated by the hyper- 
surface E* of general type and rederive the kinetic equa- 
tions for this case in Sect. II. In Sect. Ill the modification 
of the collision terms is found and the relation between 
the system obtained and the Boltzmann equation is dis- 
cussed. The correct equations of motion for the "hydro 
+ cascade" approach are analyzed in Sect. IV. 

II. DRIFT TERM FOR SEMI-FINITE DOMAIN 

Let us consider the two semi-finite domains, "in" and 
"out", separated by the hypersurface S* which in (3+1) 
dimensions will be parameterized as t = t*(x) — Xq(x). 
The distribution function 4>i n {x,p) for t < t*(x) belongs 
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to the "in" domain, whereas (j) ou t(x,p) denotes the dis- 
tribution function of the "out" domain for t > t*(x). 
Throughout this work it is assumed that the initial con- 
ditions for <fii n (x,p) are given, whereas the initial condi- 
tions for 4> ut(x,p) are not specified yet and will be the 
theme of the subsequent discussion. For simplicity we 
consider a classical gas of point-like Boltzmann particles. 

Similarly to Ref. (^] we derive the kinetic equations 
for (pi n (x,p) and <j) ut{%,p) from the requirement of par- 
ticle number conservation. Therefore, the particles leav- 
ing one domain (and crossing hypersurface £*) should be 
subtracted from the corresponding distribution function 
and added to the other one. Now we consider the closed 
hypersurface of the "in" domain, Ax 3 , which consists of 
two semi-planes o~ti and <Jt2 of constant time tl and t2, 
respectively, that are connected from tl to t2 > tl by the 
part of the boundary T,*(tl,t2). The original number of 
particles on the hypersurface o~ti is given by the standard 
expression || 



particles is conserved, i.e. N 2 = N% — N* osa + N* in , one 
can use the Gauss theorem to rewrite the obtained in- 
tegral over the closed hypersurface Ao; 3 as the integral 
over 4-volume Ax 4 surrounded by Ax 3 



J d 4 x^l p>* <9 P 4> m (x,p) 



A.i 



p v 



(f>in(x,p) - 4>out{x,p) 



e(- P v dz v ). (5) 



S*(tl,t2) 

Note that in contrast to the usual case ||, i.e. in the 
absence of boundary £*, the r.h.s of Eq. (||) does not 
vanish identically. 

The r.h.s of Eq. (||) can be transformed further to the 
4-volume integral in the following sequence of steps. First 
we express the integration element dH^ via the normal 
vector n* as follows (dx^ > 0, for j = 1, 2, 3) 

dt*(x) 



dT,^ = n* fl dx 1 dx 2 dx 



dx^ 



-(1-<W, (6) 



dS 



where d£ M is the external normal vector to o~n and, 
hence, the product p^dY,^ < is non-positive. It is clear 
that without collisions these particles can cross either hy- 
persurface <7t2 or T,*(tl,t2). The corresponding numbers 
of particles are as follows 

N 2 =JdE^p^ <f> in (x,p), (2) 
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j dH^ pf* 4> in {x,p) Q{p v dH v ) . (3) volume V# = Ax 4 U of Ax 4 and vf . Finally, with the 
J P heln of notations 



where <5 M „ denotes the Kronecker <5-function. Then, using 

tl 

identity J dt5(t — t3) = 1 for the Dirac (5-function with 
ti 

tl < t3 < t2, we rewrite the r.h.s. integral in (||) as 

dE^--- = J d 4 x 8(t - t*{x)) n* ■ ■■ , (7) 

£*(il,*2) V* 

where the 4-dimensional volume V£ is a direct prod- 
uct of the 3- and 1-dimensional volumes E*(il,t2) and 
(t2 — tl), respectively. Evidently, the Dirac ^-function 
allows us to extend integration in (R) to the unified 4- 



S*(tl,t2) 

The O-function in the loss term (^|) is very important be- 
cause it accounts for the particles leaving the "in" domain 
(see also discussion in For the space-like parts of 

the hypersurface £*(tl, t2) which are defined by negative 
sign ds 2 < of the element square, ds 2 = dt*(x) 2 — dx 2 , 
the product p v dY> l/ > is always positive and, therefore, 
particles with all possible momenta can leave the "in" 
domain through the E*(il,<2). For the time-like parts 
of E*(il,i2) (with sign ds 2 > 0) the product p h '<ffi„ can 
have either sign, and the 8-function cuts off those parti- 
cles which return to the "in" domain. 

Similarly one has to consider the particles coming to 
the "in" domain from outside. This is possible through 
the time-like parts of the hypersurface E*(tl,i2), if par- 
ticle momentum satisfies the inequality ~p v dS v > 0. In 
terms of the external normal cE M with respect to the 
"in" domain (the same as in (Q)) the number of gained 
particles 



help of notations 

e> = e(t-**(x)) ; 



= i -e- 



(8) 



it is possible to extend the l.h.s. integral in Eq. (||) from 
Ax 4 to Vy. Collecting all the above results, from Eq. (||) 
one obtains 
d 3 p 



~tl ^ 



prn 



Q(-p"nt)5(t-t*(x)). (9) 



Since volumes Ax 4 and Vy are arbitrary, one obtains the 
collisionless kinetic equation for the distribution function 
of the "in" domain 



6< p^ <pin(x,p) = 

4>in{x,p) - cf>out{x,p) 



p^n 



e(-p'X) s(t-t*(x)). (io) 



N* gain = IdS 



d 3 p 



,,-ir ^ 4>out(x,p) Q{-p v dS v ) (4) 



Similarly one can obtain the equation for the distribution 
function of the "out" domain 



E*(tl,t2) 

is, evidently, non-negative (compare it with contribution 
of feed-back particles in [(j]). Since the total number of 



©> P M 4> out (x,p) = 



p»n* 



(X,p) - (j>out(x,p) 



@{p v nl)8{t-t*{x)) (11) 
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with the same normal vector n* as in Eg. (|l0|). Note the 
asymmetry between the r.h.s. of Eqs. ( fLOf ) and (fill): for 
the space-like parts of hypersurface E* the source term 
with 0(— p"n*) vanishes identically because p v n* v > 0. 
This reflects the causal properties of the equations above: 
propagation of particles faster than light is forbidden, 
and hence no particle can (re)enter the "in" domain. 

III. COLLISION TERM FOR SEMI-FINITE 
DOMAIN 

Since in the general case 4>j n (x, p) ^ 4>out{x,p) on 
£*, the r.h.s. of Eqs. ([Toj ) and cannot vanish si- 
multaneously on this hypersurface. Therefore, functions 
= 8<|s* + and ©i, = 0>| s . ^ do not van- 
ish simultaneously on E* as well. For definiteness it is 
assumed that 

©*< = ©^ = e(o) = ~ , (12) 

but the final results are independent of this choice. 

Now the collision terms for Eqs. ([l0|) and ( pi] ) can 
be readily obtained. Adopting the usual assumptions on 
the distribution functions P), one can repeat the stan- 
dard derivation of the collision terms || and get the de- 
sired expressions. We shall not recapitulate this standard 
part, but only discuss how to modify the derivation for 
our purpose. First of all, one has to start the derivation 
in the Aa; 4 volume of the "in" domain and then extend 



Ax U Vf similarly to 



it onto the unified 4-volume Vj 4 
the preceding section. Then the first part of the collision 
term for Eq. (|Io|) reads as 



c I < (x, P ) = e 
. , i 



Ir 



{lG[4>in,4>in] ' 

D 9 P MPx) 



h [<t>ir. 



•<f>in}) , 



w, 



(13) 
(14) 



A , <t> B ] = - / D 9 P <t> A {p) <t> B {Pi) W p>PiWi , (15) 



where the invariant measure of integration is denoted as 

-* 3 — _ dV £Pi 



is the transition rate 



D 9p = iPl O^. "FL and w . , 
pU pni p/u P-P± \P P 1 

in the elementary reaction with energy-momentum con- 
servation given in the form p^ +Pj = p ,f2 +P\ ■ The r.h.s. 
of ( |l3| ) contains the square of @<-function because the 
additional 0< accounts for the fact that on the bound- 
ary hypersurface £* one has to take only one half of the 
traditional collision term (due to Eq. ( |12| ) only one half 
of E* belongs to the "in" domain). It is easy to under- 
stand that on E* the second part of the collision term 
(according to Eq. (|l2])) is defined by the collisions be- 
tween particles of "in" and "out" domains 

C< 7 (x,p) = 6<6> [Ia^inAout] - lL[<t>in,<f>out]) ■ (16) 

Combining results @, @ and @, we obtain the ki- 
netic equation for the semi-finite "in" domain 

e< p" a M <k n (x, P ) = c^p) + cH(x t p) + 



<pin(x,p) - <t>out{x,P) 



6(-pX) 5(t-t*(x)). (17) 



The corresponding equation for the "out" domain 
6> p^ <t) out (x,p) = Cl(x,p) + C"(x,p) + 



(Pin (x,p) - (f>out(x,p) 



®(p»nl)5{t-t*{x)) (18) 



can be derived similarly. In (|lj) we used the evident 
notations = Q% (Ig [pout , <j>owt] ~ h [<t>out , <?W] ) and 

Cy = ©<©> {layout, <j>in] - lL[<t>ovct,<t>iri\)- 

For the continuous distribution functions on E*, i.e. 
<A™t|s* = ^mls*, the source terms in r.h.s. of Eqs. ( |l7j ) 
and ( |l8| ) vanish and one recovers the Boltzmann equa- 
tions. With the help of the evident relations 

-d u Q<=d„, ©> =< 5(t-t*(x)) 



c 



II 



c 



II 



(19) 
(20) 



where the notation &(x, p) = 0< 4>i n (x,p)+Q> 4>out{x,p) 
is used, one can get the following result for the sum of 
Eqs. @ and @ 

p» d„ $(x,p) = -h [$,*]. (21) 
In other words, the usual Boltzmann equation follows 
from the system of Eqs. ([l7]) and (|l8|) automatically 
without any assumption about the behavior of cf> in and 
(j>out on the boundary hypersurface E* . In fact the system 
( |l7| . |l8| ) generalizes the relativistic kinetic equation to the 
case of the strong temporal and spatial inhomogeneity, 
i.e., if 4>i n (x,p) ^ 4>out{x,p) on E*. Of course, one has 
to be extremely careful while discussing the strong tem- 
poral inhomogeneity (or discontinuity on the space-like 
parts of £*) such as the so called time-like shocks pi 
because their existence may contradict to the usual as- 
sumptions adopted for distribution function. Therefore, 
in what follows we shall discuss exclusively the spatial 
inhomogeneities or discontinuities on the time-like parts 
of E* which are less restrictive because in some sense the 
equations above are delocalized in space. 

From the system (0), dl§| ) it is possible to derive the 
macroscopic equations of motion by multiplying the cor- 
responding equation with p v and integrating it over the 
invariant measure. Thus Eq. (0) generates the following 
expression (Tf = J p»p^ A (x,p)) 



©< d a ir 



d 3 p 



d3 P unii, \ 
~pP C < ( X >P) 



P 



o" P P n » 



©(- 



) S(t-t*{x)). (22) 



Similarly to the usual Boltzmann equation the momen- 
tum integral of the collision term C< vanishes due to its 
symmetries, but it can be shown that the integral of the 
second collision term C< does not vanish because it in- 
volves two different (and not identical) distribution func- 
tions. The corresponding equation for the "out" domain 
follows similarly 



©> a 



rp/MU _ 

out 



J p o VG 



(x,p) 



d 3 p 



P 



- P P U fi 



Q(p»n;)6(t-t*(x)). (23) 



Note that similar equations (with <5-like term) first were 
obtained within the relativistic hydrodynamics in pi. 
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IV. DISCUSSION 

It is clear that Eqs. @, @, (§|) and (||) remain 
valid for the finite domains as well. With their help we 
are ready to analyze the "hydro+cascade" models. In the 
TLS model the cut-off formula relates (j>i n (= hydro, Eq. 
(f22|)) and <fi out (= cascade, Eq. (|l8|)) on E* as follows 



TLS 



<d(p p n* p ) (f) out 



&(p P n p ) 0, 



(24) 




i.e., for the space- like parts of hypersurface E* these func- 
tions are identical, whereas for the time-like parts of E* 
there are no returning particles to the "in" domain. In 
this case the source term in cascade Eq. ( |l8| ) is zero, 
while the source term in hydro Eq. (|2^) does not vanish 
on the time-like parts of the boundary E*. Therefore, 
the main defect of the TLS model is not even the energy- 
momentum non-conservation, but the incorrect hydrody- 
namic equations. The absence of the 5-like source term 
breaks the conservation laws (evidently, the system 
231) obeys the conservation laws), but its inclusion 
into consideration will inevitably change the hydrody- 
namic solution of Ref . Q . The full analysis of the possi- 
ble solutions of the systems ( |i"7[|i"6| ) and ( p2| , |l8| ) requires 
a special consideration. We only mention that from the 
negative sign of the TLS source term in the r.h.s. of (|22] ) 
for equal indices v — fi one immediately can deduce that 
such a correction to the hydro equations should increase 
the degree of the fluid rarefaction in comparison with the 
standard hydrodynamic expansion. It is, therefore, quite 
possible that such a source term will generate a discon- 
tinuity between "in" and "out" domains. In the ther- 
modynamically normal media |TT| the rarefaction shocks 
are mechanically unstable. However, it is well know that 
on the phase transition boundary between QGP and HG 
the properties of the mixed phase are thermodynamically 
anomalous [pi]] and the usual rarefaction shocks are pos- 
sible. Another possibility is the occurrence of the new 
type of the discontinuity, the freeze-out shock suggested 
in Refs. where the post freeze-out state is described 
by the cut-off distribution and, hence, is very similar to 
the TLS ansatz. It is clear that in both cases the ad- 
ditional rarefaction will reduce the mean transverse size 
and the life-time of the hadronizing QGP. 

Let us consider briefly the BD approach. Since in the 
BD model the hydro and cascade distributions on E* are 
equal (f> ou t\s* — 4>i n \?.*, the corresponding source terms 
vanish in all equations. Therefore, at first glance the 
BD approach correctly conjugates the hydro and cascade 
solutions on the arbitrary hypersurface. For the over- 
simplified kinetics considered above it is so. However, 
the real situation differs essentially from our consider- 
ation. Thus, the hydro part in both the BD and TLS 
models is assumed to be in the local thermodynamic 
equilibrium, whereas the matter in the cascade domain 
can be far from equilibrium |t],|| (this was, actually, the 
main reason why both groups decided to use the cas- 
cade). Consequently, the BD transport equations for all 



hadrons are homogeneous in the hydro domain, whereas 
for the cascade domain they are inhomogcneous. Since 
the initial BD cascade distribution on the time-like parts 
of E* contains the particles returning to the fluid domain 
(j> v n* v < 0), then these particles will move towards the 
space-like parts of the hypersurface E* which are located 
inside of the light cone originated at each point of the 
time- like part of E*. Then for each hadron the inho- 
mogeneous BD cascade equation will generate a different 
distribution function than the one already obtained from 
the hydro equations (or homogeneous transport ones) on 
these space-like parts of E*. Thus, one arrives at a causal 
paradox, the recoil problem jl| , and at a mathematical in- 
consistency. 

Evidently, the inclusion of the viscosity into the hy- 
dro equations (apart from its tremendous complexity for 
the mixture of about hundred hadrons) will not solve the 
problem because for the small deviations from equilib- 
rium (and, hence, a small viscosity effect) the influence 
of the returning particles may remain essential. If, on the 
contrary, the matter in the hydro domain is far from equi- 
librium, then the usage of the hydro equations becomes 
problematic. Therefore, a more realistic way is to find 
the boundary conditions for cj)i n and <j) out on the separat- 
ing hypersurface E* form the system of kinetic equations 
(O.Eq), and then to apply these boundary conditions to 
the system ( |f7| , ^2] ) which ensures the correct treatment 
of the relativistic nuclear collisions within the frame of 
the "hydro+cascade" model. 
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